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In this paper the KTeV collaboration reports the analysis of 112.1 x lO'^ candidate Kl — > 7r+7r~7 decays including a 
background of 671±41 events with the objective of determining the photon production mechanisms intrinsic to the de- 
cay process. These decays have been analyzed to extract the relative contributions of the GP violating bremsstrahlung 
process and the GP conserving Ml and GP violating El direct photon emission processes. The Ml direct photon 
emission amplitude and its associated vector form factor parameterized as j^MiKl + (m^-m^)+2M e ) have been 

measured to be |gMi| = 1-198 ± 0.035(stat) ± 0.086(syst) and 01/03 = -0.738 ± 0.007(stat) ± 0.018(syst) GeVVc^ 
respectively. An upper limit for the GP violating El direct emission amplitude \gEi\ < 0.21 (90%GL) has been found. 
The overall ratio of direct photon emission (DE) to total photon emission including the bremsstrahlung process (IB) 



has been determined to be DE/(DE +IB) 



0.689 ± 0.021 for E^ > 20 MeV. 
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The study of the direct photon emission in K^ -^ 
7r+7r~7 decays gives insight into both the structure of 
the kaon and the sources of GP violation in this mode. 
This decay proceeds via two main processes [1,2]. The 
first of these is the inner bremsstrahlung process (IB) 
in which one of the charged pions from a GP violating 



Kr 



T'^T 



decay emits an El electric dipole photon by 



bremsstrahlung. The second process is the emission of a 
GP violating El electric dipole photon or a GP conserving 
Ml magnetic dipole photon together with the tt+tt^ pair 
directly from the primary decay vertex. The photons pro- 
duced by the IB process have a typical bremsstrahlung 
spectrum with E.^ in the Kl center of mass peaking to- 
ward zero and falling off like 1/E^, while the direct pho- 
ton emission produces an energy spectrum peaked toward 
larger E^. The Ml amplitude is expected to require a 
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perturbation theory [4] in order to incorporate the ef- 
fects of the structure of the K^ on photons emitted at the 
primary decay vertex (as opposed to the photons of the 
bremsstrahlung process emitted from the charged pions). 
In this form factor M"^ and M]^ are the mass squared of 
the p (770 MeV/c^) and K (497 MeV/c^) mesons.. 

The KTeV collaboration previously reported a mea- 
surement [3] using 8669 candidate K^ — > 7r+7r~7 decays 
accumulated in the 1996 KTeV E832 run at Fermi Na- 
tional Accelerator Laboratory which indicated clearly the 
presence of the Ml process and the need for the associ- 
ated form factor. In addition, the presence of Ml photon 
emission and the need for a form factor have also been 
demonstrated by the KTeV E799 experiment [5,6] and 
NA48 [7] measurement of the Kl -^ 7r+7r~e+e~ mode. 
While this mode differs from the Kl — * 7r+7r~7 since 
the photon is virtual converting internally to a e^e~ 



Dalitz pair in the Kl -^ 7r+7r~e"'"e~ decay, both modes 
have the same ampUtudes contributing except for the 
presence of an extra "charge radius" amplitude in the 
Kl -^ 7r"''7r~e+e~ decay. Thus we expect the same gMi 
amphtude and associated form factor to be present in the 
Kl -^ TT+TT^e+e^ decay. This has been demonstrated by 
the measurements detailed in Refs [5-7]. 

Values of |(?mi|, its form factor, and the ratio 
DE/(DE+IB) presented in this paper were determined 
by using the much larger, complete KTeV E832 1997 
Kl — > 7r"'"7r~7 data set containing 112.1 x 10'^ can- 
didate Kl -^ 7r+7r~7 decays. We have analyzed the 
Kl -^ 7r"''7r~7 decay mode using the double differential 
decay rate 
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from the model of Rcf. [8]. In this expression, 9 is the 
angle of the photon with respect to the Tr"*" in the tt+tt" 



center of mass system, 
Kl rest frame, and 



Ej is the photon energy in the 
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\Mldirect\ = |5Ml|(l + 
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(M2 - M^) + 2MkE. 
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where 5o(s = Mj^) and 5i{s = M^^) are the isospin = 
0,1 strong interaction tt+tt" phase shifts evaluated at the 
kaon mass and at the particular tt+tt" mass of a given 

Kl -^ 7r+7r~7 decay. \ri^ |e**+~ is the amplitude for the 

CP violating Kl — > tt+tt" decay. 

Note that there is no interference term between the 
El and Ml amplitudes. However, there can still be 
an interference term in the differential decay rate be- 
tween the EIbr and Eldirect amplitudes. The interfer- 
ence will generate a contribution to the E^ energy spec- 
trum intermediate in energy between the lower energy 
bremsstrahlung photons and the higher energy Ml pho- 
tons. 

The Kl — * tt^tt^j signal of 111.4 x 10'^ events above 
a background of 671± 41 events, obtained after the anal- 
ysis cuts described below, is shown in Fig. 1. Details of 
the detector and beam can be found in Ref. [10] so we 
only give a brief overview here. A proton beam with 
a typical intensity of 3 x 10^^ delivered in a 20 sec- 
ond spill every minute was incident at an angle of 4.8 
mr on a BeO target producing two nearly parallel Kl 
beams, one of which intercepted a Kg regenerator and 
the other of which remained a "vacuum" beam. The 
data for the Kl -^ n^n^j measurement were obtained 
from the "vacuum" beam decays. The configuration of 
the KTeV E832 vacuum beam and detector consisted of 



a vacuum decay tube, a magnetic spectrometer with four 
drift chambers, photon vetoes, a Cesium Iodide (Csl) 
electromagnetic calorimeter, and a muon detector. 

Approximately 4.3 x 10* events were extracted from 
two track triggers [10] by requiring that the two tracks 
to pass track quality cuts and form a vertex with a good 
vertex x^- These tracks were also required to have op- 
posite charges and E/p < 0.85, where E was the energy 
deposited by the track in the Csl, and p was the momen- 
tum obtained from magnetic deflection. Showers chosen 
as photons are required to be far from pion showers and 
to have a transverse shower shape consistent with elec- 
tromagnetic showers. Only photons with E^ > 20 MeV 
in the tt^tt~j rest frame were included in this analysis. 

To reduce backgrounds arising from other types of Kl 
decays in which decay products have been missed, the 
candidate 7r"*'7r~7's were required to have transverse mo- 
mentum Pf relative to the direction of the Kl be less 
than 2.5 x 10^^ GeV^/c^ and M^+^+^, the invariant mass 
of the 7r+7r~7 system, to be 490 MeV/c^ < M„^^ < 
506 MeV/c2. 
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FIG. 1. TT 7r^7 invariant mass for events passing all 
Kl -^ TT^TT^j physics cuts except for the M^+^-^ cut. 
Crosses are data and the solid line is the fit to the background 
components 

The major background to the Kl — > 7r+7r~7 mode 
was due to an accidental calorimeter cluster in coinci- 
dence with a Kl — > it fi^v decay in which the muon 
was misidentified as a pion. This background was sup- 
pressed by the muon detector identification as well as 
the Pi and Mt^.^^ cuts. A smaller background was due 
to Kl —>■ TT^e^v decays in which there was an acciden- 
tal photon, and the electron was misidentified as a pion. 
This background was suppressed by electron E/p identi- 
fication and Pj? and M^^t^^ cuts. The Af^+^-^ spectrum 
shapes due to the K^z and K^'^ backgrounds were similar. 



A still smaller background to the Kl -^ 7r+7r~7 mode 
was Kl — > tt+tt^tt'^ in which one of the photons from 
the 7r° decay was not detected in the Csl calorimeter or 
the photon vetos. To reduce the Kl — > 7r+7r~7r'' back- 
ground, the longitudinal momentum (P|)^o of all can- 
didate Kl —f 7r+7r^7 events was calculated (under the 
assumption that the events were really K^ —>■ tt'^tt^'tt'^) 
in the frame where the tt+tt" momentum was transverse 
to the Kl direction. In this frame, (P^)^o is > 0.0 (except 
for resolution effects) for Kl -^ 7r"'"7r~7r° decays. In con- 
trast, the Kl -^ 7r+7r~7 decays should have {PDttO < 0. 
The requirement -0.10 < (P|)^o < -0.0055 GeVVc^ 
together with the P^ and Mj^^r-y niass cut, suppressed 
the Kl — > tt'^tt^tt'^ background. 

Hyperon decays such as A — > p7r~ plus an accidental 
photon with the proton misidentified as a n'^ , or S -^ 
Att^ with a misidentified proton and one of the tt" pho- 
tons missed, were determined to contribute a few events. 
Other sources of background such as Kl -^ tt+tt^ coin- 
cident with an accidental photon or Ks -^ 7r+7r^7 pro- 
duced in the neutral beam production target were com- 
pletely negligible. 

The magnitude of the remnant background after all 
cuts was determined by a fit (see Fig. 1) of the sideband 
regions above and below the Kl mass peak to shapes ob- 
tained from a Monte Carlo of the backgrounds leading to 
an estimated total background of 671 ± 41 events. The 
best estimate of the composition of this background is 9% 
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and the remainder due to the other minor backgrounds 
mentioned above. The 112.1 x 10'^ candidate events, in- 
cluding the estimated 671 events of background, were an- 
alyzed in a likelihood fit based on equations 1 and 2. The 
likelihood was a function of the two independent vari- 
ables 9 and E^, the values of the fit parameters ai/a2, 
\gMi\ and \gEi\ and nominal values from the PDG [11] 

for the other model parameters such as ry^ The strong 

interaction phase shifts of the tt+tt^ system are taken 
from Ref. [12]. The likelihood was calculated using a 
Monte Carlo event sample generated with nominal val- 
ues of the fit parameters, traced through the spectrome- 
ter undergoing multiple scattering, bremsstrahlung, and 
secondary decays. The resulting events are then recon- 
structed using the same reconstruction code as was used 
on data. These reconstructed Monte Carlo events are 
then reweighted with a new set of fit parameters using 
the Kl — * 7r+7r~7 matrix element of Ref. [8] and a likeli- 
hood is calculate for the new parameters. The maximum 
likelihood fits to the two independent variables cos,6 and 
-E-y are shown in Fig. 2a) and Fig. 2b) respectively. 

Possible systematic uncertainties in ai/a2, |5a/i| and 
\gEi\ due to disagreements between data and Monte 
Carlo simulations were investigated by varying analysis 
cuts and observing variations in these fit parameters. In 
addition, the momentum spectrum of the 7r"'"7r~7 system 
observed in the Kl — >• tt+tt" decays has been adjusted 



to agree with the Kl momentum spectrum observed in 
K^ -^ 7r+7r~ mode and the data has been refit after the 
adjustment. Any differences between ai/a2, \gMi\ and 
\gEi\ before and after the final adjustments were taken 
to be a systematic error due to uncertainty in the kaon 
beam momentum spectrum. Systematics due to uncer- 
tainties of parameters such as rj-^ , and the strong in- 
teraction phase shifts 5o,i that were not determined by 
the fit were studied by varying each parameter over ±lcr 
of their published values and observing the variation of 
0-1 / 0,2-, I^Mij, and |(7bi|. Bremsstrahlung radiation from 
the pions was studied using the PHOTOS program [13]. 
This radiation could lead to 7r^7r~77 final states in which 
one of the photons is not observed causing shifts of the 
kinematics of the original 7r7r7 decay. Possible system- 
atic effects due to the non-orthogonality of drift chamber 
planes were also studied. Final overall systematic errors 
in ai/a2, |5a/i|, and \gEi\ were obtained by adding the 
individual errors in quadrature. Table I lists the non-zero 
systematic uncertainties of ai/a2 and \gMi\, and |5bi|. 
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FIG. 2. Likelihood fit to the two independent variables 
in the Kl rest frame: a) the angle 9 between the tt"*" and 
the 7 in the tt'^tt~ center of mass and b) the photon energy 
spectrum E^ in the Kl rest frame. The components of the 
photon energy spectrum and the cos^ spectrum due to the 
bremsstrahlung and Ml direct emission processes are shown. 

The results, including systematic errors of the mea- 
surement of the Ml direct photon emission amplitude 



and the attendant vector form factor, are ai/a2 — 
(-0.738±0.007(stat)±0.018(syst)) GeVVc^ and |5Afi| = 
1.198±0.035(stat)±0.086(syst). These measurements are 
in good agreement with the measurements of Ref. [3,5-7] 
(see Fig. 3). After incorporating the systematic errors, 
an upper Hmit of \gE\\ < 0.21 (90% CL) was obtained. 

Using the result for \gM\ \ and its associated form factor 
and taking \gE\\ to be equal to zero, the ratio of direct 
to total photon emission in K^ —* n'^n^j decay was 
determined by integrating the Ml and bremsstrahlung 
processes over 9 and Ej (for E~^ > 20 MeV) to be 



DE/(DE+IB)= 0.689 ± 0.021. 
agreement with Ref. [3]. 



This result is in good 



TT+TT^e+e" decays [6]. 
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FIG. 3. 90% CL contours of qmi vs. axjai for vari- 
ous experimental measurements; 90% CL results from the 
Kl -^ 'K^Tr~'y mode from this paper (5-shaded contour) with 
the 68% CL contour also shown (white contour); For compar- 
ison, we show the results from the Kl -^ iT'^-K~e^e~ mode 
for NA48 data (3-dotted contour) of Ref. [7], for KTeV 1997 
data (2-dashed contour) in Ref. [5], and for KTeV 1997-^1999 
data (4-solid contour) of Ref. [6], and for axja-i from earlier 
KTeV 96 Kl -^ 7r"'"7r~7 data (1-Iight gray vertical region) of 
Ref. [3]. 

In conclusion, this paper presents the best mea- 
surements achieved to date for the Ml direct pho- 
ton emission form factor parameters \gM\\ = 1.198 ± 
0.035(stat)±0.086(syst)andai/a2-0.738±0.007(stat)± 
0.018(syst) GeVVc^ in the Kl -^ tt+tt--/ and Kl -^ 
7T^TT~e'^e~ decay modes. These measurements are in 
good agreement with our previous measurement of ai/a2 
using the 1996 KTeV Kl -^ tt+tt-j data [3] and with 
our measurements of I^miI and axja-i using the 1997 
and 1999 KTeV Kl -> Tr+Tr-e+e" data [5,6] and with 
NA48 results [7] from Kl — > 7r+7r~e"'"e~. We have also 
determined an upper limit \gE\\ < 0.21 (90%CL) for 
CP violating El direct photon emission in the Kl -^ 
7r+7r^7 mode consistent with that measured using Kl -^ 



[9] 



[10] 

[11] 
[12] 

[13] 



Rev. D46, 



D48, 



L.M. Sehgal and M. Wanninger, Phys 

1035(1992); ibid. D46, 5209(E)(1992). 

P. Heiliger and L.M. Sehgal, Phys. Rev 

4146(1993). 

A. Avati-Harati et aZ.,Phys. Rev. Lett. 86, 761(2001). 

Y.C.R. Lin and G. Valencia, Phys. Rev. D37, 143(1988) 

A. Avati-Harati et al, Phys. Rev. Lett. 84, 408(2000). 

E. Abouzaid et al, hep-ex/05080010. 

A. Lai et al. Eur. Phys. Jou. C30, 33(2003). 

L.M. Sehgal and J. van Leusen, Phys. Rev. Lett. 83 

4933(1999). 

J.K. Elwood, M.B. Wise, and M.J. Savage, 

Rev. D52, 5095(1995); J.K. Elwood et al. 

D53, 2855(E)(1996); J.K. Elwood et al, ibid. _ 

4078(1996). 

A. Avati-Harati et al, Phys. Rev. D67, 012005(2003): 

ibid. D70, 079904(2004). 

S. Eidelman et al, Phys. Lett. B592, 010001(2004). 

S. Pislak et at, Phys. Rev. Lett. 87, 221801(2001). 

E. Barberio and Z. Was, Comput. Phys. Commun. 79 

291(1994). 



Phys. 
ibid., 
D53, 



Source 


\9mi\ 


ai/a-z 


\9El\ 


Differing initial MC parameters 


0.0093 


0.0021 


0.013 


Kaon Beam Momentum Uncertainty 


0.0031 


0.0004 


0.005 


Background uncertainty 


0.0355 


0.0067 


0.045 


Pion bremsstrahlung 


0.0326 


0.0140 


0.097 


Non-Orthogonality of chambers 


0.0402 


0.0013 


0.009 


Physics cut variations 


0.0463 


0.0056 


- 


Fitting resolution 


0.014 


0.0056 


0.024 


Ej,cos6 resolution 


0.023 


0.0042 


0.038 


77+ _ uncertainty 


0.0171 


0.0014 


- 


So phase uncertainty 


0.0111 


0.0021 


- 


5i phase uncertainty 


0.0053 


- 


- 


Total Systematic Error 


0.086 


0.018 


0.117 



TABLE L Contributions to the systematic errors (67% CL) 
for ai/a2, \gMi\, and \gEi\. 



